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Shape-dependent thermodynamics and non-local hydro- 
dynamics are argued to occur in dissipative steady states of 
driven diffusive systems. These predictions are confirmed by 
numerical simulations. Unlike power-law correlations, these 
phenomena cannot be explained by a hypothesis of "critical- 
ity". Instead, they require the effective Hamiltonian of the 
system to contain very long-range potentials, making the in- 
variant probability measures formally "non-Gibbsian" . 



themselves. They are expected to occur in physical drift- 
diffusion systems, e.g. electrolytes and semiconductors. 

Correlations of power-law type are predicted in such 
systems by linear fluctuating hydrodynamics In a 

DDS with one species of (unit-charged) particle at con- 
stant mean density p, in contact with a heat bath at tem- 
perature T and in an applied electric field E, the equation 
in Fourier space for density fluctuations Sp(k, t) is 



PACS Numbers: 05.70.Ln, 66.30.Hs, 82.20.Mj, 64.60.Lx Sp{k, t) = [ic E -k - k-D E -k]5p(k, t) - ik-S3(k, t). (1) 



Power-law decay of correlations are generic in dissi- 
pative steady-states of open, driven systems with con- 
servation laws In equilibrium systems such power 
laws can be due to two different mechanisms: the inter- 
action potentials in the system Hamiltonian can them- 
selves be long-range power-laws (e.g. dipolar) or else 
the potentials can be short-ranged but the system may 
be at a critical point ||. The latter circumstance has 
prompted the view that dissipative, noncquilibrium sys- 
tems are attracted without any tuning of parameters to 
a critical state, or exhibit so-called "self-organized crit- 
icality" (SOC) [Q. We shall show here that such an 
interpretation, taken literally, is not true in an impor- 
tant class of such systems. As we shall explain, an effec- 
tive Hamiltonian may be introduced to characterize the 
noncquilibrium statistics. We then exhibit phenomena 
that can be explained only by the presence there of ef- 
fective power-law potentials of very long range. Such an 
explanation for power-law correlations was earlier pro- 
posed in driven diffusive systems Q and, more recently, 
in a shaken granular flow Because the underlying 
dynamics are local in these systems, this mechanism 
may be justly termed "self-organized long-range interac- 
tions" (SOLPJ). We shall demonstrate in the class of sys- 
tems considered two closely related phenomena: shape- 
dependent thermodynamics and non-local hydrodynamics. 
These phenomena could not occur if the effective Hamil- 
tonian were short-ranged but critical. In fact, the in- 
duced potentials must have such an extreme long-range, 
many-body character that the noncquilibrium measures 
are formally non-Gibbsian. (For an excellent review of 
the relevant notions, see 0.) The indicators of the 
non-Gibbsian nature, shape-dependent thermodynamics 
and non-local hydrodynamics, have great importance in 



Here Ce(p) = %rG°) is the drift velocity, a density deriva- 
tive of the conduction current. The latter is given, at 
low density p and small field-strength E, by Ohm's law, 
Je = er-E = /O/x-E, with <r the conductivity tensor and 
/j, the mobility tensor. De(/)) is the diffusion tensor in 
the conducting state, while S3(k, t) is the Fourier compo- 
nent of the current noise with zero mean and covariancc 
(SJ(k,t)53(k',t')) = 2L E (p)<5 3 (k + k')5(t - t'). In the 
limit E — y 0, Le(p) coincides with the Onsager matrix 
L = a-ksT. It is easy to solve Eq.([j]) for the static struc- 
ture function SE(k) = limt_>oo(<$p(k, t)Sp(— k, t)}: 



5 E (k) = 



k-L E -k 
k-D E -k 



(2) 



See Inverse Fourier transforming gives a long-range 
density-density correlation S'e(i') ~ r~ d in dimension d, 
when Le and De are not proportional. This will gener- 
ally be true for E ^ 0. Such power-law correlations have 
been verified in numerical simulations of simple driven 
lattice gas (DLG) models Correlations of this type 

are expected more generally in locally driven/damped 
conservative systems without detailed balance. For ex- 
ample, r~ d correlations of similiar origin have recently 
been predicted and verified computationally in the ho- 
mogeneous cooling state of rapid granular flows . Our 
results have implications for all such systems. 

We illustrate our points with the DLG models P,^2[. 
These were originally introduced as models of solid elec- 
trolytes, or superionic conductors [ p3||l4| ]. However, un- 
like physical drift-diffusion systems of charged particles, 
these models have only short-ranged dynamical interac- 
tions and are thus perfect for our theoretical objectives. 
We shall comment later how the results carry over to re- 
alistic DDS with Coulombic interactions. The particles 
of the DLG model live on a cubic lattice Z d . Assuming 
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hard-core exclusion, the occupancy i] Xyt e {0, 1} for each 
site x £ Z d and time t > 0. The evolution of the config- 
uration is via a Kawasaki exchange dynamics, specified 
by the rate ce (x, y; 77) for exchange of occupancy of near- 
est neighbor sites x, y in the configuration 77, i.e. for the 
transition 77 — > r) xy . In addition to assuming that rates 
are functions of occupancies at sites within a finite range 
of {x, y}, the main assumption is local detailed balance: 

c E (x, y; 77) = c E (x, y; n xy ) exp [-/3 {H{rf y ) - H(rj) 

+E.(x-y)(r ?x -7, y ))], (3) 

for some short-ranged lattice-gas Hamiltonian if (77), e.g. 
an Ising model. The condition (^) encourages particles to 
hop in the direction of the electric field E and, in infinite 
volume, sets up an irreversible steady-state with a mean 
current. In fact, these models have space-ergodic, ho- 
mogeneous, and time-invariant measures (J- p ,e,0 for each 
density p £ [0, 1], expected to be unique at small /?. The 
question arises whether these measures in infinite volume 
are "Gibbsian" for an effective Hamiltonian 



(4) 



AcZ d 



with some set of many-body potentials $,4 depending 
on spins rj x at sites x £ A C Z d . If it exists, this will 
generally not be the same as the short-ranged Hamilto- 
nian H(ri) used in defining the dynamics. It turns out 
that almost any reasonable measure (with local densi- 
ties) is "Gibbsian" if one permits extremely long-ranged, 
many-body potentials and, indeed, there are a myriad 
of physically inequivalent such Hamiltonians! (E.g. see 
fL5| |, Theorem V.2.2(a)). To guarantee uniqueness and 
other standard properties of usual Gibbs measures, the 
condition of absolute summability is required: 



A3x 



< OO 



(5) 



for all x € Z d , where ||$aJJoo = sup^ ^^(r})). Following 
rather common practice 0, we shall agree here to call 
only measures with the latter property "Gibbsian" . It 
is a rigorous theorem of Asselah jl6| (or Appendix B of 
H ) that the following alternative holds: either all space- 
ergodic, invariant measures of the DLG are Gibbsian with 
absolutely summable potential or else none of them are. 
It may also be proved that in the domain of analyticity, 
no absolutely summable power-law potential can produce 
a correlation cx r~ d (Theorem 1 of ||), such as is ob- 
served in the DLG. Thus, under the first alternative, one 
is led to conclude that the Gibbs measure must be critical 
to account for the observed correlation decay. Alterna- 
tive #1 for our purposes may thus be termed SOC. On 
the other hand, in alternative #2 the potentials are non- 
summable (= long-ranged), so that this case corresponds 
to the SOLRI scenario. To support the latter, we make 
some key comparisons with long-ranged systems. 



An analogy was already remarked a few years ago 
between the DDS and dipolar systems ||. Of course, 
dipole spin-spin correlations are also oc r~ d even at 
high-temperature. This is consistent with our point of 
view, because the dipole potential just misses being abso- 
lutely summable (and is thus "non- Gibbsian" according 
to our criterion!) An important consequence of this non- 
summability was early recognized (Tt]], namely, that the 
thermodynamics of dipolar systems is shape- dependent. 
This situation arises because the dipole potential energy 
sums are only conditionally convergent and hence may 
lead to different values depending upon the order of sum- 
mation, i.e. the shape, at least at nonzero field H|. In 
uniformly-magnetized (=high field), ellipsoidally-shaped 
samples of dipolar materials the shape dependence of 
thermodynamic free-energy functions is simply param- 
eterized |fl9| , in good agreement with experiment po| . 

It was argued in ^ that a similiar shape-dependence 
occurs in DDS. The thermodynamic functions of interest 
are the "pressure" £>e and the "Hclmholtz free-energy" 
/e- The former is defined by the thermodynamic limit 



Pe(p, 0\p*) 
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lim 

A^z d /3\A\ 



(6) 



where A is a sequence of lattice volumes converging to 
Z d , Na(ti) is the number of particles within A for the 
configuration 77, and the average (-)p^,/3,E is with respect 
to the invariant measure /i Pi ,/3,E of the DLG for refer- 
ence density /?*. The "Helmholtz free energy" /e is 
then introduced by the Legendre transform /e(/0, (3\p*) = 
supJfip ~ pe(p, /3\p*)]- We include the reminder 
of the reference density and the E to indicate the strength 
of the applied electric field. Since the measures here 
are for irreversible steady-states with Ohmic dissipa- 
tion, these are not usual free energies. They coincide 
with the equilibrium free-energies in the limit E — > 0. 
The physical interpretation of Je{p, P\p*) is as an "ex- 
cess dissipation function," i.e. as the total energy dis- 
sipated per volume by an external field to change the 
density to p from its reference value p* (in addition to 
the Ohmic dissipation intrinsic to the reference state) 
||. The argument for shape-dependence is that the 
"susceptibility" (essentially, the isothermal compressibil- 
ity) may be written both in terms of the free-energy, 



Xe(p,P) = P 



and also in terms of the 



structure- function, via the limit xe = limt^o SE(k). 
However, the latter limit is indeterminate when the struc- 
ture function has the form in Eq.(||) and depends upon 
the wavenumber vector direction k along which the limit 
is taken. Thus, the free-energy itself must be shape- 
dependent, by the same argument as for dipole systems. 

To test this prediction we have performed a Monte 
Carlo simulation of the DLG on a periodic square S x S 
lattice with Ising Hamiltonian H = _ ^( x y )My 
where (x, y) denotes nearest neighbor sites, for which the 



2 



(inverse) critical temperature is (3 C « 0.31 To stay 

well within the single phase (high temperature) regime, 
we used [3 — 0.2, E — 10.0 and reference density p» = 0.5. 
We have determined the thermodynamic functions for 
rectangular subblocks A of the S x S system, in which 
various aspect ratios of the sides of the rectangles were 
chosen. The pressure was evaluated by a double limit. 
First the infinite volume limit was obtained by a linear 
extrapolation in 1/S — > on the Monte Carlo average 
{ e P» N *) s in the steady-state with S = 64, 128,256,512. 
The thermodynamic limit in Eq.(^) for the pressure was 
then evaluated by a second linear extrapolation in the 
inverse volume 1/|A| of the subblock going to zero. The 
largest subblock edge in this second extrapolation had 
length 18. The Legendre transform to the free-energy was 
then carried out. The results are shown in Figure 1 for 
subblocks with aspect ratios of 3:1 and 1:3 for sidelcngths 
parallel and perpendicular to the field, respectively. Er- 
ror bars reflect both statistical deviations in independent 
runs and the double extrapolation procedure. The two 
functions are clearly distinct. We see that the DLG, con- 
sidered as a model of a current-carrying electrochemical 
cell, has well-defined free-energies but the results depend 
upon the shape of the cell! 

This shape-dependence leads, however, to a rather seri- 
ous puzzle about the hydrodynamic behavior of the DDS. 
The general problem is to describe how an initial smooth 
density profile relaxes to a constant density in the driven 
steady-state. In Q a nonlinear hydrodynamic equation 
was derived by the formal method of noncquilibrium dis- 
tributions to describe this irreversible process. Density 
fields varying on a length-scale of the order of e _1 com- 
pared with the lattice distance were formally shown to 
evolve by a drift- diffusion equation, 



p(r,t) = -V- 



%(p)-/?l e (p).v(^) 



(7) 



over times of order e~ 2 . This equation has the "On- 
sager form", with j E (p) the conduction current, L E (p) 
the Onsager coefficient matrix, and T E [p] the free-energy 
functional. Explicit analytical formulae were given in 
H for each of these quantities, e.g. a Green-Kubo for- 
mula for the Onsager matrix. These formulae are ex- 
act even at high field strengths E, although they may 
be difficult to evaluate concretely. The free-energy func- 
tional in [f| was nominally given by the local expression 
■7"e[p] = / dr/ E (p(r)). However, as observed there, such 
a form is indeterminate. Since the free-energy depends 
upon the limiting shape, which value is to be used? 

We can now resolve this issue. The free-energy func- 
tional actually shown in Q to be relevant to hydro- 
dynamics is given by a Legendre transform T E [p] = 
J drp(r)p(r) — "P E [p] of the pressure functional 

V E [p] = \ lim e d log(exp[V/3 M (ex)? 7x ]) p ^ 3 , E . (8) 

' X 



Simple computation then gives g^jTS 
5 2 V E 



p* and 



Sp(r)5fi(r') 



= P I™ £ d [(V[e-^r]V[e-^r'])p,,E - pi 

= [3S E {v-v'). 



(9) 



It follows that V E [p] = f drp*(i(r) + § f dr Jdr'S E (r - 
r')/i(r)/i(r') + 0(p 3 ), and the Legendre transform yields 

Te[p] = ip~J dr J ^'S^(T-v')Sp(r)Sp(T') + 0(Sp 3 ) 

(10) 

where 5p(r) = p(r) — p* and S^ 1 is the operator in- 
verse of S E . In other words, 5'g 1 (k) = k-D^-k/k-L^-k. 
Hence, the inverse kernel S^(r) is a r for large r, 
too. We see that the hydrodynamic equation of the DDS 
at finite field strengths E thus must have an explicit, 
severely non-local form. [H. Spohn has emphasized to us 
that Eq.(^|) linearized about the homogeneous state of 
density p* will still be local if Eq.(|^) holds, since then 
L E (p*):VV / dr'Vtr-rOMr') = D E (p*):VVSp(r).] 
However, with the free-energy functional in ( |ITj| ) replac- 
ing the local expression, all results of M remain valid: 
the H-theorem, the fluctuation-dissipation theorem, etc. 

Such nonlocal hydrodynamic behavior should also be 
present in equilibrium systems with long-ranged interac- 
tions, e.g. dipolar hard sphere systems or ferrofluids. In 
one such system, the Kawasaki lattice gas with a long- 
ranged Kac pair-potential, there is a rigorous result |2lf| 
that a hydrodynamic equation of the same form as (|7j) 
holds, with a similar nonlocal expression for the free- 
energy as in (|l0|), simply replacing S , £ 1 (r) by the kernel 
J(r) of the Kac potential. We see again a very striking 
and fruitful analogy between dissipative, driven systems 
and equilibrium systems with long-range interactions. 

Our results verify that the SOLRI scenario holds in 
our model, and not SOC. For Gibbsian measures with 
summable potentials there is no shape-dependence of 
thermodynamics, such as we observe here, even at the 
critical point ( |15[ |, Theorem 1.2.5). Boundary condi- 
tions play a role below the critical point in the phase co- 
existence region in determining which of multiple phases 
will occur, but even then the free energies are indepen- 
dent of the phase. The ordinary thermodynamic limit, 
with no shape-dependence, remains valid directly at the 
critical point. Likewise, the dynamics of short-ranged 
systems in the coexistence region is expected to be de- 
scribed by the Cahn-Hilliard dynamics, which is local. It 
is actually a little perplexing how to interpret the SOC 
point of view that nonequilibrium steady-states are al- 
ways "critical" , when these are observed themselves to 
undergo continuous phase-transitions at sharp values of 
temperature and/or density. Such transitions occur both 
in the DLG |12| and in granular flow 0], not to mention 
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dipole systems p2] . Away from the transition point both 
DDS and dipole systems have a finite correlation length £ 
which characterizes the crossover from a critical power- 
law oc r~' d ~ 2+? ') at intermediate range r <C £ into the 
asymptotic power-law oc r~ d at long-range r » £ Jl2| , p2[ . 

It is most interesting to consider the implications of our 
model calculation for real systems. We expect that the 
main conclusions concerning shape-dependence and non- 
locality will hold for physical drift-diffusion systems, such 
as fluid or solid electrolytes, semiconductors, and, at a 
more mesoscopic level, colloidal suspensions [ p3| . In real 
charged-particle systems there is the added complication 
of a dynamics which is itself long-ranged, via Coulombic 
interactions. However, these are expected to be Debye- 
screened and effectively short-ranged. The drift-diffusion 
equations were justified long ago for nonequilibrium pro- 
cesses in electrolyte solutions at low density and small 
fields within Debye-Huckel theory [24j] . The equations 
are of the same form as those we have considered, sim- 
ply generalized to multiple ionic species. The effects 
of the Coulomb interactions are calculable and can all 
be incorporated into an effective Onsager matrix, with 
cross-species terms due to ionic-cloud distortion and elec- 
trophoresis. The effects considered in our work corre- 
spond to contributions to the invariant measures and 
thermodynamic potentials in at least the E 2 power of 
the field strength. It would be interesting to make a 
theoretical estimate of the order of magnitude of such ef- 
fects in electrolyte solutions. Perhaps the most accessible 
predictions are the long-ranged correlations themselves, 
which could be observed in light-scattering experiments 
similiar to those carried out on simple fluids subject to a 
temperature gradient (see pj| for a recent review) . Also 
of possible practical interest are the implications of a non- 
local hydrodynamics for granular flow. 
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FIGURE CAPTION 

Figure (1.) Free energy as a function of density for 
E = 10, (3 = 0.2, p» = 0.5 for aspect ratio 3:1 (triangles) 
and 1:3 (squares). 
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